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Abstract 

^ , Simulation studies are undertaken for the system Ag/Au(100) by means 

'^ ■ of grand canonical Monte Carlo applied to a large lattice system. The 

interactions are calculated using the embedded atom model. The forma- 
tion of adsorbed Ag phases of low dimensionality on Ag(lOO) is investi- 
gated and the influence of surface defects on the shape of the adsorption 
isotherms is studied. The results of the simulations are discussed in the 
jrt ■ light of experimental information available from electrochemical measure- 

ments. 
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O. 1 Introduction 

New horizons have opened concerning the understanding of metal electrodeposi- 
^ ' tion processes in the last decade. This has been achieved to a large extent due to 

^D , the use of electrochemical techniques like cyclic voltammetry along with in-situ 

nanoscopy techniques. Among these, scanning tunnelling microscopy (STM) 
has allowed the observation in direct space and real time of the formation and 
growth of a new phase during the electrocrystallization of a metal from an elec- 
trolytic solution. In this context we can mention the so-called underpotential 
deposition (upd), which implies the deposition of a metal M on the surface of 

(^ , an electrode or substrate S of different nature at potentials more positive than 

those predicted from the Nernst equation. A recent and extensive discussion on 
this phenomenon can be found in a book of Budevski, Staikov and Lorenz [T], 

^\^ • and the theoretical aspects have been considered in reviews by one of us [H H] . 

H , One way to quantify upd is through the so-called underpotential shift 

A0[/P£), which is related to the difference of the chemical potential of M ad- 
sorbed on S* at a coverage degree 0, say ijl{Mq/ S), and the chemical potential 
of M in the bulk phase, say fi{M/M), through the following equation: 

A^UPD = —[fi{M/M) - ti{Me/S)] (1) 
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where z is the charge of the ion M in the solution and cq is the elemental 
charge. Criteria for the experimental estimation of ^(^upd have been given in 
the literature by Kolb et al.[4] and Trasatti [5]. Concerning the modelling of 
upd as a thermodynamic process, different interesting approaches have been 
taken by several authors in the 70's [TJ [SI [3 [HI H] ■ The most widespread idea is 
to envisage this phenomenon as related to the formation of different adsorbed 
phases, involving in some cases, phases of low dimensionality that differ in their 
properties and can be detected as separate entities by a given experimental tech- 
niques, i.e. peaks in a voltammogram. The interpretation of upd as a phase 
transition has been employed by Blum et al [5] to develop a model for the un- 
derpotential deposition of Cu on Ag(lll) in the presence of sulfate anions and 
has also been recently brought to a more general context by Lorcnz and cowork- 
ers [in]. According to these ideas, the dimensionality of the phases involved in 
the upd phenomena is between and 2. The OD phases correspond to the 
adsorption of isolated metal atoms, the adsorption on kink sites, the adsorption 
on vacancies of the first lattice plane of the substrate or on some other type 
of point defects. The ID phases are related to the adsorption at the border of 
steps or monoatomic terraces and the 2D phases are related to the formation of 
monolayers and submonolayers of adsorbate building compressed, (1x1) or even 
expanded structures. 

Figure 1 illustrates the most common situations experimentally observed, 
which will also be considered in the present simulations studies. 

It can be observed that relative abundance of OD and ID phases is deter- 
mined by the presence of defects on the electrode surface. Its existence has been 
related to the existence of small shoulders in the voltammograms, as postulated 
by Siegentaler et al. [TT] and Lorenz and coworkers [H]. In the case of the 
2D phases that are formed on the flat terraces or defect-free surfaces there is 
clear experimental evidence for the existence of compressed [13j . compact and 
expanded structures [Il[in]- The formation and growth of a compact, defect- 
free bi-dimensional structure can be described as a first-order phase transition 
that takes place at the electrode/solution interface and its occurrence should be 
characterized by a discontinuity in the adsorption isotherm, with a concomitant 
very sharp peak (more properly speaking a Dirac delta function) in the voltam- 
metric profiles f^. Concerning this point, it is remarkable that although upd 
has been often related to the presence of first-order phase transitions, the real 
existence of these phenomena in experimental systems is still not completely 
free of controversy, since discontinuities in the experimental isotherms are usu- 
ally not obtained when specific anion adsorption is avoided [T] . We arrive thus 
at the somewhat paradoxical situation where experimentalists usually perform 
their experiments under conditions where the clean adsorption process of metal 
adsorption is actually coupled to the even more complex process of anion coad- 
sorption. Higher order phase transitions should be characterized by the presence 
of rather wide and flat voltammetric peaks [9J , though there is no experimental 
evidence for the occurrence of this type of phenomena in electrochemical sys- 
tems T. In this respect, computer simulations could also deliver enlightening 
information. Furthermore, computer simulations play a fundamental role for the 
comprehension of surface phenomena ,21|. Many approaches can be taken for 
the study of surface-related problems. Among others, we can mention molecular 



dynamics simulations, Monte Carlo simulations based on continuum Hamilto- 
nians (ofF-lattice) and techniques in which the surface is modelled through a 
lattice Hamiltonian (on-lattice simulations) |22j . 

The present work deals with the computer simulation of Ag deposition of 
Au(lOO) by means of the Grand Canonical Monte Carlo method and the inter- 
atomic potential given by the embedded atom method (EAM), which has been 
previously employed by some of us to study electrochemical systems [Tl]. We 
developed a simulation scheme that combines the simplicity of the lattice model 
of the surface with the many body properties of the metallic binding. Adsorp- 
tion isotherms under different temperature conditions were obtained in order 
to analyze the existence of a first-order phase transition, and the influence of 
different surface defects on the behaviour of the isotherms is also investigated. 

2 Some experimental facts 

Cyclic voltammetry and STM in-situ experiments performed on the systems 
Ag/Au(hkl) have shown the presence of different adsorbed phases, depending 
on the applied electrode potential [TUl [13 [13 [2U] . In the particular case of 
Ag adsorption on Au(lOO), the experiments have revealed the existence of an 
expanded, quasi hexagonal c(-\/2 x 5-\/2)it!45°Ag phase at high overpotentials. 
The existence of this expanded structures can hardly be explained in terms of 
energetic considerations, so that the entropic contributions to the free energy 
of the system should be considered. Furthermore, the presence of specifically 
adsorbed anions may in turn play an important role in stabilizing this type of 
entities [T5j , so we shall leave the consideration of this type of structures for fu- 
ture work. On the other hand, at intermediate underpotentials the decoration 
of monoatomic steps corresponding to the formation of a one-dimensional phase 
was observed, while at underpotentials close to Ag bulk deposition, two dimen- 
sional islands were built on terraces and defect-free zones of the substrate. STM 
measurements with lateral atomic resolution showed that these islands present a 
quadratic symmetry, suggesting that they may be considered as a (1 x 1) struc- 
ture commensurate with the Au(lOO) substrate. Finally, it must be added that 
long time polarizations showed some indication for the beginning of a surface 
alloying process whose mechanism has so far not been explained. 

3 Model and simulation method 

3.1 Lattice model 

Lattice models for computer simulations are of widespread use in studies of nu- 
cleation and growth, because they allow dealing with a large number of particles 
at a relatively low computational cost. In principle, it must be kept in mind that 
continuum Hamiltonians should be much more realistic in those cases where epi- 
taxial growth of an adsorbate leads to incommensurate adsorbed phases j23j or 
to adsorbates with large coincidence cells. On the other hand, the use of fixed 
rigid lattices restricts enormously the number of possible configurations for the 
adsorbate and its use may be justified on the basis of experimental evidence or 
continuum computer simulations that predict a proper fixed lattice geometry. 
In the present case, we have strong evidence from this latter type of simulations 



within the canonical Monte Carlo method 24] that indicates that at least one of 
the phases present during Ag underpotential deposition on Au(lOO) possesses a 
pseudoniorphic structure. In fact, our continuum MC simulations showed that 
a Ag monolayer adsorbed on Au(lOO) spontaneously acquired a (1 x 1) coin- 
cidence cell in agreement with the experimental finding at low underpotentials 
already mentioned in section 2. For this reason, we shall employ here a lattice 
model to represent the square (100) surface lattice in a Grand Canonical Monte 
Carlo simulation. Besides Ag adatoms, Au adatoms are allowed on the surface 
at different coverage degrees, building a variety of surface structures with the 
purpose of emulating some of the most common surface defects. Linear sweeps 
of the chemical potential are performed to obtain the adsorption isotherms and 
study the influence of the surface defects on their shapes. First of all, we neglect 
the effect of the presence of solvent molecules. This approximation should not 
be critical as long as the partial charge on the adatoms is small, thus minimiz- 
ing the ion-dipole interactions. Second, we also neglect all kinds of anion effects 
that may coadsorb during the metal deposition process. In the case of other 
systems, the adatom-anion interactions have been shown to be very important 
[17j , playing a decisive role in determining the energetics of the system. This 
has been recently analyzed by some of us in thermodynamic terms [H] . In the 
case of the present system, however, thin layer twin electrode experiments gave 
no indication for any change in the extent of anion adsorption upon building 
the adsorbed monolayer. 

Square lattices of different sizes with periodical boundary conditions are 
used in the present work to represent the surface of the electrode. Each lattice 
node represents an adsorption site for a Ag or a Au atom. The former may 
adsorb, desorb or jump between neighboring sites, while the latter may only 
move on the surface like the Ag atoms do. In this way, our model corresponds 
to an open system for one of its components, that is, Ag. This has physical 
correspondence with the setup of the electrochemical experiment, where only 
the metal in equilibrium with its cations in solution may dissolve or be deposited 
in the potential range considered. 

Concerning the Au adatoms, some considerations must be made regarding 
the existence of surface defects. These atoms may in principle move freely on the 
surface and minimize the free energy of the system by a number of mechanisms. 
For example, isolated Au atoms may heal defects through their incorporation to 
a defective cluster, or small Au islands may dissolve to join large ones, like shown 
in previous simulation work by Stimming and Schmickler |25j . However, as long 
as these processes are slow enough as compared with the processes related to 
the Ag adatoms (and this is found to be the case), the whole process can be 
envisaged as the adsorption/desorption behaviour of Ag on a Au(lOO) surface 
with defects. Thus, different Au structures can be imposed as initial conditions 
for each simulation. In the present case, they range from islands of different 
sizes and shapes obtained by means of simulated annealing techniques, up to 
completely random distributions. All this was undertaken in order to emulate 
some of the defects that can be found on a real single crystal surface, like kink 
sites, vacancies, isolated substrate atoms, steps, etc. 



3.2 Interatomic potential 

A very important feature to be taken into account when comparing the results 
of a simulation with experiment is the quality of the interatomic potentials used 
to perform the simulations. Several methods have been developed to calculate 
the total energy of a many-particles metallic system, with a computational effort 
comparable to that of a pair potential [26'. Among these models, denominated 
broken-bonds or bond-cutting methods it is worth mentioning the embedded 
atom method(EAM) ^7\, the N-body potentials of Finnis and Sinclair [IS], the 
second-moment approximation or Tight-binding (TB) |26j and the glue model 
(GM) [29]. In this work we use the embedded atom method [27] because it is 
able to reproduce important characteristics of the metallic binding that cannot 
be obtained using simple pair potentials. 

The EAM considers that the total energy Utot of an arrangement of N par- 
ticles may be calculated as the sum of energies Ui corresponding to individual 
particles 

N 

Utot =yu, (2) 



i=l 



where Ui is given by 



C/» = F.(pm) + ^E^'^('^^^-) (3) 
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Fi is denominated embedding function and represents the energy necessary to 
embed atom i in the electronic density ph,i at the site at which this atom 
is located. ph,i is calculated as the superposition of the individual electronic 
densities Pi{rij) : 

PKt = yPi{rij) (4) 

Thus, the attractive contribution to the EAM potential is given by the em- 
bedding energy, which accounts for many-body effects. On the other hand, 
the repulsion between ion cores is represented through a pair potential Vij{rij), 
which only depends on the distance between the cores r^ : 

Tij 

Zi{rij) may be envisaged as a sort of effective charge, dependent on the nature 
of the particle i. The EAM has been parametrized to fit experimental data like 
elastic constants, dissolution enthalpies of binary alloys, bulk lattice constants 
and sublimation heaths [57] . Pair functionals have been widely used for surface 
diffusion studies and adsorption of metals on metallic surfaces [TH [301 [H] ■ 

3.3 Grand Canonical Monte Carlo 

One of the most appealing characteristics of Grand Canonical Monte Carlo 
(pVT/MC) is that the chemical potential p, is one of the independent variables, 
as it is the case in many experimental situations. In the electrochemical case, 
we can mention low-sweep rate voltammetry, where the electrode potential can 
be used to control the chemical potential of species at the metal/solution inter- 
face. This technique offers a straightforward way of obtaining the adsorption 



isotherms provided the sweep rate is low enough to ensure equihbrium for the 
particular system considered. Another interesting aspect of fjiVT/MC is that 
a better and faster equilibration of the system is possible due to local density 
changes originated by fluctuations in the number of particles. 

Our 2D system is characterized by a square lattice with M adsorption sites. 
We shall label each adsorption site with a number 0, 1 or 2, depending on 
whether it is empty or occupied by a gold or a silver atom respectively. Thus, 
the configuration of the system will be defined through the components ni, 
"■2,' ■ •,'T-Af of an Af -dimensional vector Tt and its partition function may be 
written as: 

M M 

E{ni,fi2,M, T) = ^ exp/3(^i ^ 5i,„^ + /^2 JI ^2,„J exp(-/3f/(7i^)) (6) 

where fii and /i2 are the chemical potentials of gold and silver at the interface, 
<5i,ni and ^2, Hi are Kronecker delta functions acting on the components of the 
vector ~rt and U{~rt) is the energy of the system in state Tt . The probability of 
obtaining a given configuration will be given by the corresponding Boltzmann 
factor: 



p / , ,r ^N exp/?(^i ^fio Ji,„^ + At2 E^=o '52,n.) exp(-/?[/(Tf )) 
P^(Mi,A^2,A/,T)- -———- (7) 

Following the procedure proposed by Metropolis and coworkers [32], the 
acceptance probability for a transition from state ft to it is defined as: 

W^Tf^Tf' =^i^(l.^) (8) 

n 

so that detailed balance is granted. 

In our ixVT / MC simulation we shall allow for three types of events : 

1. Adsorption of a Ag atom on a lattice site selected randomly. According 
to eqnsl7]and[Sl the acceptance probability W^^f^^^f' for this event is: 

Wjf^^^ = min(l, exp(/3(Mi - AC/))) (9) 

where AC/ = C/(Tf' ) - C/(Tf ). 

2. Desorption of a Ag atom selected randomly. This probability is 

Wj^^j^, = min(l, exp(/3(-/ii - AC/))) (10) 

3. Motion of a Ag or a Au atom from the lattice site where it is adsorbed to 
one of its four nearest neighbors. The latter is selected randomly and the 
acceptance probability for this event is equivalent to that employed in the 
canonical ensemble, since the number of particles remains constant. We 
have: 

Wj^^jf. - min(l, exp(-/3AC/)) (11) 



Even when in a grand canonical simulation, events of type 3 (motion or dif- 
fusion of lattice particles) are not strictly necessary, their presence is justified by 
the fact that a smaller number of Monte Carlo steps(MCS) must be employed for 
the equilibration of the system[32]. The calculated equilibrium thermodynamic 
properties do remain unaltered by this choice. 

In order to satisfy detailed balance, the a priori probability for the creation 
of a Ag atom, ac^ must be equal to the probability of its anhilation ad and 
independent of its motion probability am- 

The algorithm devised for an elementary step of the simulation is the fol- 
lowing: 

1. Selection of the elementary process (adsorption, desorption of diffusion) to 
occur. This is made by means of the generation of a random number e 
uniformly distributed in the interval [0, 1]. 

2. If < e < {ac + ad), a random site is chosen on the lattice. 

• If the occupation number of the site is 0, an attempt is made to create 
a Ag atom at that site, according to the acceptance given by eqn. [SI 

• If the occupation number of the site is 2, an attempt is made to 
destroy the atom adsorbed at that site, according to the acceptance 
given by eqn. 1101 

3. If {ac + ad) < e < 1 , an attempt is made to move each of the particles of 
the system according to the following algorithm: 

• One of the nearest sites to the particle is selected. If the occupation 
number is 0, a movement to the new position is attempted according 
to eqn. [TTJ Otherwise, the next particle in the list is considered. 

4. Back to point 1, thus completing one MCS. 

Within this procedure, the relevant thermodynamic properties are then ob- 
tained as average values of instantaneous magnitudes stored along a simulation 
run. For example, the average coverage degree of the silver atoms (O)^ at a 
given chemical potential fii after a simulation with MCSe equilibration steps 
and a total number of MCS steps will be given by: 

MCS 

(e(M:))., = ,,^^_,,^^^ ^^£^9(,,)... (12) 

where the instantaneous value of silver coverage degree, 0(/ii)Ag,i , is defined 
from the number of Ag atoms NAg,i and the number of Au atoms Nau present 
at the interface at the time step i 

e(^^), ^^^^d^^ (13) 

Since our main goal is to obtain adsorption isotherms for several Ag coverage 
degrees and different surface structures of Au, we chose a simulation method 
that does not deliver information on the kinetics of nucleation and growth. In 
this respect, the Kinetic Monte Carlo(KMC)[351[35 technique should be apphed 
and will be the subject of future work. 



3.3.1 Algorithm employed for the calculation of energy differences 

In a MC simulation the determining step concerning the speed of the calculation 
is the evaluation of the total energy of the system U{~rt), which is determined 
by the type of interatomic potential employed. 

According to the algorithm described in the previous section, any new state 
of the system is accepted with a probability that involves the Boltzmann fac- 
tor of the energy difference between the initial and the final state, that is 
AC/ = U{lt') — U{~rt). One of the main advantages of the lattice model is 
its simplicity, since it fixes the distances between the adsorption nodes, thus 
reducing the energy values that the system can take to a discrete set. Further- 
more, the potentials used are short ranged, so that a very important simplifying 
assumption can be made for obtaining A[/. The point is to consider the ad- 
sorption(desorption) of a particle at a node immersed in a certain environment 
surrounding it, as shown in Figure 2 . The adsorption site for the particle is 
located in the central box, and the calculation of the interactions is limited to 
a circle of radius R. Then, the adsorption energy for all the possible configura- 
tions of the environment of the central atom can be calculated previous to the 
simulation. The S adsorption sites within the circle are labelled arbitrarily with 
numbers between 1 and 5* (Figure 2) and their labels along with the occupation 
number may be employed to identify the configuration in terms of a number 
/ expressed on the basis of 3. Generally speaking, for a configuration char- 
acterized by the occupation numbers ni , n2 , . . . ng , the corresponding value 
of / is calculated as / == X]i=i ^j 3'^^*^^ . For example, for an environment 
with 8 neighbors {S — 9) such as that shown in Figure 2a, the configuration 
Til = 2, n2 = 2, jij, = 1, n^ — 1, n<^ = 0, hq — 0, n^j — 0, n^ — is indexed 
according to our convention with / = 6156. With this method all the adsorp- 
tion/desorption energies of Ag and Au on Au(lOO) are tabulated, so that during 
the MC simulation the most expensive numerical operations are reduced to the 
reconstruction of the number / that characterizes the configuration surrounding 
the particle on the adsorption node. Computationally speaking, / is nothing 
but the index of the array in which the energy is stored. 

In the present work, the adsorption/desorption energies were calculated for 
environments with R equal to the distance between second and third nearest 
neighbors, as shown in Figure 2 , using as a substrate a five-atomic planes thick 
gold slab. The results for the energy are shown as a histogram in Figure 3 , 
where the frequency of appearance of a given energy is shown as a function of 
the adsorption energy. This is a very important information, since this figure 
contains the distribution of energy changes that an atom may undergo during 
the simulation. A more detailed discussion on the information that can be 
extracted from the histograms will be given in section 4.1. 

We must finally mention a further assumption in our calculations. In all the 
structures generated the adsorbates were located at a vertical distance from the 
first substrate plane equivalent to the distance between neighboring Au(lOO) 
lattice planes. This should not be a very important approximation, since the 
difference between the Au and Ag lattice constants amounts only 0.01 A. 



3.3.2 Simulated annealing 

Simulated annealing techniques have often been used to obtain minimal energy 
structures or to solve ergodicity problems. A suitable way to implement them 
is through the canonical Monte Carlo method at different temperatures. The 
simulation is started at a very high initial temperature To, of the order of 10^ 
K, and the system is later cooled down following a logarithmic law: 

Tf = ToK^"^"'" (14) 

where T/ is the final temperature, Ncydes is the number of cooling steps and 
ii' is a positive constant lower than one. A few hundreds of thousands of MCS 
are run at each temperature in order to allow an extensive exploration of the 
configuration space and the simulation stops when Tf is reached. Although 
experimental annealing treatments (like flame annealing) are used to get clean 
and structurally well defined surfaces, in the present case we have used the com- 
putational homologue to obtain the more stable surface defects of Au at 300 K. 
In the case of experiments, this procedure has been sometimes questioned and 
we will see here that the type of structures obtained are indeed very sensitive 
to the cooling rate employed. 

4 Results and discussion 

4.1 Some considerations on the adsorption/desorption en- 
ergies of Ag 

The information provided by the histograms shown in Figure 3 is very valuable 
since they contain all the possible adsorption/desorption energies that a Ag atom 
may undergo during the simulation. For the purpose of the present analysis, we 
roughly distinguish between three types of configurations: those in which the 
central atom (box 13 in fig. 2b) is surrounded only by other Ag atoms(or empty 
sites), those in which the central atom is only surrounded by Au atoms and 
those in which the environment may contain Au or Ag atoms indistinctly. The 
two former families involve 4096 configurations each, while the latter consists of 
531441 configurations, including the two previous ones. In can be observed that 
in the two former cases (continuous and dotted histograms respectively) the 
energy distribution is multimodal, they being qualitatively very similar to each 
other. However, the histogram corresponding to Ag adsorption in the presence 
of Au atoms is flatter and is shifted towards more negative adsorption energies. 
This suggests a certain affinity of Ag for a Au environment. From this point 
of view, the adsorption of Ag in the presence of Au surface defects should be 
favoured with respect to the adsorption on a perfect surface. This speculation 
will be corroborated by the MC simulations. 

The pulse-plot in Figure 3 corresponds to the distribution of the adsorption 
energies in the mixed environment for the 12-neighbors system. The energy 
distribution looks still multimodal but considerably smoother than the previous 
ones. Thus, in general terms it can be concluded that the presence of gold atoms 
surrounding the adsorbed Ag atoms as surface defects will widen and enrich the 
energy spectrum of the adsorbate. 



4.2 Ag deposition on a perfect Ag(lOO) single crystal sur- 
face. 

Our initial work was devoted to the deposition of Ag on a defect-free Au(lOO) 
surface, performing iJ,VT / MC simulations at 300 K. The systems consisted in 
(40 X 40) and (100 x 100) lattices, using the adsorption/desorption energy tables 
calculated as described in Section 3.3.1. 

Positive and negative chemical potential sweeps were made in the range 
between — 2.7ey and —3.2eV, involving simulations with a length of 10® to 
10^ MCS , depending on the system size. In order to quantify underpotential 
deposition, the deposition of Ag on Ag(lOO) was also simulated. The chemical 
potential for bulk Ag deposition /^5i„ /^„ was determined from these runs. 

In all cases the adsorption isotherms show a very abrupt change of the cov- 
erage degree. Fig. 4a, at a characteristic chemical potential /i^ ,^^ , that mainly 
depends on the cutoff radius employed for the calculation of the energy and to 
a lower extent on the system size. The /^j4„/^„ values along with other rele- 
vant information are reported in Table 1. The drastic change observed in the 
adsorption isotherms points towards the existence of a first-order phase transi- 
tion in which a two-dimensional phase (an adsorbed monolayer) is built through 
a nucleation and growth process. According to the values of chemical poten- 
tials reported in Table 1, this occurs at lower /i values than those expected 
for bulk Ag deposition (mJ4„m„ = —2.865 eV), leading to an underpotential 
shift A(j)irpD — iJ-Ag/Ag ^ fJ'Ag/Au that Can be estimated to be between 0.06 eV 
and 0.11 eV . The isotherm corresponding to Ag adsorption on Ag(lOO) is also 
shown in Figure 4a. 

Concerning the chemical potential sweep, it was performed in two different 
fashions. In the first, for each new ^ value characterizing the simulation, an 
initial value of coverage degree of Ag, Oaq, was set, that was 0, 0.5 or 1. This 
delivered three types of isotherms that corresponded to the different initial con- 
ditions. In the second type of simulations, the final configuration at a given 
jj, was employed as starting point for the next simulation. The first type of 
simulations was employed to study hysteresis effects, that in our system were 
found to be important at 300 K (Figure 4b). As shown below, they are related 
to the absence of surface defects since it is precisely on these defects where the 
formation of the new phase starts. The simulations undertaken with the second 
scheme showed essentially the same characteristics as those with the first. 

The effect of temperature was also investigated and it was found that at high 
temperatures the sharp behaviour of the Q — /i isotherms changes to a smooth, 
langmuirian-like behavior, as can be seen in Figure 5. From these curves we can 
infer that at 300K the present system is well below its critical temperature. We 
can also estimate this critical temperature using the equation of the 2-D Ising 
model: 



2fcln(V2- 1 



where k is the Botzmann constant and w is the interaction energy between 
nearest neighbours. We have estimated a value of w by inserting a Ag adatom 
in different environments. The values obtained were in the range —0.27 eV < 
w < —0.18 eV so that the critical temperature is estimated to be in the range 
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1200 K < Tc < 1800 K. These values should only provide a rough estimation, 
since the Ising model does not consider the many body effects present here. 

4.3 Influence of surface defects on the adsorption isotherms 

A wide variety of defects may exist on a single crystal surface, to an extent that it 
depends among other factors on the temperature, the chemical environment and 
the method by which it was generated. We therefore performed an important 
number of simulations to elucidate the effect of the most common surface defects 
on the Ag adsorption isotherm. The starting configuration for the Au defects on 
the surface was obtained from a simulated annealing calculation and a 100 x 100 
lattice and R ~ 3 were used in all cases. The length of the simulations ranged 
between 10^ and 10^ MCS, with an adequate number of equilibration steps. 

4.3.1 Generation of surface inhomogeneities by means of simulated 
annealing 

It has been demonstrated experimentally that after a flame annealing treat- 
ment the Au(lOO) surface is reconstructed to give a quasi-hexagonal structure 
[35j . This reconstruction is lifted upon a potential sweep at a potential value 
that depends on the nature of the electrolyte involved. Since the hexagonal 
structure is more compact than the (1 x 1) structure, a number of islands are 
built on the electrode surface that accommodate the remaining atoms. In the 
case of the experimental studies of Ag deposition on Au(lOO), STM and AFM 
in-situ images indicate that previous to Ag deposition the substrate consists in 
terraces separated by monoatomic steps, while the corresponding images with 
lateral atomic resolution show a quadratic structure characteristic for an unre- 
constructed surface pUfT^ . 

Since the atomic excess expected after the lifting of the reconstruction is 
about 13 % of the atomic density of the Au(lOO) single crystal surface, we 
generated Au islands that cover the single crystal surface with a coverage degree 
Oau — 0.1. Thus, simulated annealing was performed with 1.000 Au atoms 
distributed on a 100 x 100 lattice at an initial temperature Tq of 10.000 K. The 
system was later cooled down to a final temperature T/=300 K following the 
cooling rate given by eqn. [H] with K= 0.8. This produced 16 cooling steps 
between the initial and the final temperature. In order to produce surfaces with 
different densities of defects, the number of Monte Carlo steps allowed at each 
temperature was varied, obtaining the different final structures shown in Figure 
6 . It can be noticed that when the number of MCS at each temperature 
is increased, the islands turn more compact, with well defined borders. This 
is equivalent to a low concentration of defects such as kinks and monoatomic 
steps. On the other hand, if the cooling is faster, the number of islands per unit 
surface increases and their borders turn more irregular. Thus, the lower the 
number of MCS allowed at each temperature, the farther the final state will be 
from equilibration. Table 2 shows the density of kink sites pk and the density 
of sites located at the border of the monoatomic steps ps at different number 
of MCS. The pk was calculated as -j^ where M^ is the number of kink sites 
and M ~ N'^ is the total number of lattice sites(lO.OOO). Analogously, ps was 
calculated as -^ where Ms is the number of sites located at the border of the 
monoatomic steps. The relationship perimeter/area of the islands, calculated 
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from the number of Au atoms located at the border and the number of atoms 
located inside the islands, is also given in the table. All the values refer to the 
final configuration of the system. 

4.3.2 Deposition on kink sites 

The first stage in the deposition process is given by the adsorption of Ag atoms 
on point defects. This can be envisaged as the formation of a 0-dimensional 
phase that occurs at chemical potentials more negative than the chemical po- 
tential for the formation of the Ag monolayer. Although the most commonly 
observed point defects on the surface of a monocrystalline electrode at 300 K 
are kink sites, their density at this temperature may be affected by different 
factors, like the presence of a condensed phase in contact with the metallic sur- 
face [T] . In the present work we considered that the initial number of kink sites 
and the corresponding roughness of the monoatomic terraces of the substrate 
only depends on the cooling rate employed in the simulated annealing process. 
The initial structure of the stepped Au surface was set according to the results 
obtained in section 4.3.1 and sweeps of chemical potential were undertaken be- 
tween —3.30 and — 3.02ey . Figure 7a shows isotherms obtained with 5 of 
the 16 configurations resulting from the simulated annealing process. The gold 
structures considered correspond to the final state resulting from simulations 
with 20, 320, 5120, 81920 and 655360 MCS and range from polyhedron shaped 
islands with different roughness up to a more or less random distribution of very 
small irregular islands. A coverage degree due to adsorption on kink sites O^ 
may be defined considering that they form a sublattice of M^ adsorption nodes 

n'i 
at a given stage of the simulation. According to this, we defined 0fc = "m'j 

where A'^ is the number of Ag atoms adsorbed on kink sites. It can be clearly 
observed that none of the isotherms present an abrupt change in the coverage 
degree, showing a rather langmuirian shape that indicates that cooperative ef- 
fects are absent in this type of adsorption. Another remarkable feature that 
is relevant for the discussion below is the fact that for increasing density of is- 
lands and the concomitant increasing roughness of their borders, the adsorption 
isotherms are shifted towards more negative chemical potentials. 

It is clear that the kink sites subject of the present analysis are not all 
equivalent. Furthermore, many of them are actually not independent since the 
distance between them is lower than the cutoff radius for the interactions and 
they can even move in the course of the simulation due to the mobility allowed 
to the Au atoms. It may thus be interesting to compare the present isotherms 
for kink sites with others where the equivalence and independence of the ad- 
sorption nodes is granted a priori by considering an ideal ensemble of identical 
and independent sites. With this purpose we built such lattices with the most 
commonly found kink sites, resulting from our simulated annealing procedure. 
The most commonly found types of kink sites are shown in figure 8 . With 
the adsorption energies corresponding to each of these sites we constructed the 
corresponding Langmuir isotherms for the coverage degree Q^ ""^ as a function 
of the chemical potential /i. These isotherms were then used to compose what 
we shall denominate an " ideal " adsorption isotherm for Ag adsorption on kink 
sites Q^kink/ 1^ according to the following. Let us denote with p\, pi and p| 
the densities of each of the types of kink sites calculated from the simulated 
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annealing results aceording to: 



Ml 



with Ml denoting the number of kink sites of type i and Mk = M^ + M| + M|. 
From the three isotherms Qi ""^//^, ©2 ""^/m ^-^d O3 ""^/a*, we defined the "ideal 
" adsorption isotherms for adsorption on kink sites as: 

©S(m) - pI^i"'"' + pI^^"^' + p|e3^""^ (15) 

Ideal kink sites isotherms for the five imperfect structures selected are shown 
in Fig. 7b. Comparison with Fig. 7a indicates that the more compact the islands 
the closer the behaviour of the real isotherms to the ideal ones. On the other 
hand, for increasing roughness of the island borders, the real isotherms shift 
their inflection point towards lower chemical potentials. This larger deviation 
from ideality in the case of the rougher structures is understandable, since in 
this case the distance between kink sites diminishes and they no longer behave 
as independent adsorption sites. 

Table 3 summarizes the chemical potentials at the inflection point (fimfi) 
for real and ideal isotherms. Comparison with the value for the deposition of 
a Ag monolayer on a perfect surface shows that the presence of point defects 
may shift the threshold for underpotential deposition in the negative direction 
as much as 0.2 V. 

It is also interesting to compare these values with that corresponding to the 
filling of Au surface monovacancies with Ag atoms. In this case we obtained 
H'^nfi — ^3.53 eV. However, these defects are expected to have a too low 
concentration at the working temperatures in electrochemistry so their relevance 
may be considered as negligible for practical purposes. 

4.3.3 Deposition on the border of monoatomic steps 

The second stage of the present electrodcposition process involves the formation 
of one-dimensional phases through the decoration of monoatomic steps. In all 
cases this process was found in the ^ range between —3.05 eV and ~-2.95eV, 
that is, after deposition on point defects and before the formation of the 2-D 
phase. The corresponding simulations were run with the five systems selected 
above, that is, a 100 x 100 lattice partially covered with Au islands and the 
chemical potential was run between —3.32 eV and —2.90 eV. Figure 9 shows the 
adsorption isotherms for the surface configurations selected. Analogously to the 
case of point defects, we defined here a coverage degree related to border sites, 
Qstep, referred to the total number of adsorption sites located at the border of a 
monoatomic step in a given configuration of the simulation. If these isotherms 
are compared with those of figure 7a, it can be noticed that the roughness of the 
Au islands plays the opposite role concerning the stability of the Ag deposit. 
That is, the increasing roughness of the islands shifts the isotherms towards more 
positive chemical potentials. This fact may be in principle explained in terms 
of the considerations made in section 4.1 about the influence of Au atoms on 
the energy spectrum of Ag adsorption. Rougher islands provide adsorption sites 
with a lower coordination of the Ag border adatom with the Au atoms of the 
island, so that the adsorption energy of the adsorbates is in the average shifted 
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towards more positive values. However, this is a rather simpHfied explanation 
and a more detailed study in terms of the surface configurations should be 
undertaken. 

Finally, it is remarkable the fact that as the compactness and size of the 
islands increase, the isotherms become steeper, indicating some sort of cooper- 
ative effect in the formation of the 1-D phase. 

4.3.4 Formation of Ag monolayers in the presence of surface defects 

The formation of the 2-D phase in the presence of Au islands was finally in- 
vestigated in connection with the other processes described above. Chemical 
potential sweeps were undertaken in the range [— S.SOeF < ^ < — 2.90ey], 
with the observation of adsorption of Ag atoms on kink sites [—3.30 eV < /i < 
—3.02 eV], decoration of monoatomic steps, [—3.05 eV < ^i < —2.95 eV] and 
the final formation of the monolayer[ —2.95 eV < /i < —2.90 eV]. Figure 
10 shows the final state at different /i for one of the five systems considered 
previously(MCS' = 5120). 

If we compare the isotherms for this systeni(Fig. 11) with those correspond- 
ing to adsorption on perfect surfaces (Fig 4a), it is clear that the presence of 
surface defects introduces additional stages in the deposition process due to 
the formation of 0-D and 1-D phases. It is remarkable the smoothing of the 
isotherms in the neighborhood of the deposition potential of the monolayer, this 
being enhanced for rougher substrates. It is also observed that as the number 
of islands per unit surface increases, the chemical potential for the formation of 
the monolayer is shifted negatively. In all cases the formation of the 2-D phase 
is built in a coalescence regime starting from the Au islands. This can be noted 
in the snapshots of Figure 10. 

All the deposition stages described occur at chemical potentials negative with 
respect to Ag bulk deposition {fi'X„/An ~ —2.865 eV), thus within the under- 
potential deposition range. The underpotential shifts Acj^upu for the different 
processes thus range between 0.4 eV and 0.03 eV. 

5 Conclusions 

The deposition of Ag on Au(lOO) was studied in the absence and in the presence 
of surface defects. These were generated by means of a simulated annealing 
procedure using different numbers of Monte Carlos steps [MCS) per cooling 
cycle. It was found that a lower number of {MCS) generates a higher number of 
Au islands per unit surface, which are characterized by a higher perimeter/area 
rate. On the other hand, a lower cooling rate generates more compact, bigger 
Au islands. 

Ag deposition always takes place at chemical potentials lower than bulk de- 
position, thus leading to a positive underpotential shift. This is also supported 
by first-principles calculations performed by one of us [36]. The adsorption 
isotherms obtained indicate that Ag adsorption should occur in three well de- 
fined stages: first, adsorption on kink sites, second, adsorption at the border of 
monoatomic steps and finally formation of the monolayer or 2-D phase. In all 
cases the presence of inhomogeneities shifts the adsorption isotherms towards 
more negative chemical potentials thus increasing the underpotential shift with 
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respect to that predicted for a perfect surface. Furthermore, the isotherms are 
smoothed in the neighbourhood of the critical chemical potential, this effect 
being more pronounced for stronger disordered surfaces. 

The results of our simulations support to a large extent the experimental 
results obtained for this system, fundamentally concerning the existence of a 
positive underpotential shift and the strong affinity of Ag for Au evidenced in 
the trend to build surface alloys. More detailed experimental studies concerning 
the influence of surface defects on the voltammograms and adsorption isotherms 
would be very helpful to check the present calculations and to improve the 
modelling of the system Ag/Au(100). One point to be considered is that of 
finite size effects of the simulation cell, especially in the case of adsorption on 
kink sites for imperfect surfaces obtained at a low cooling rate. 
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7 Tables 



R \N X N 


40 X 40 


100 X 100 


2 


-2.975 


-2.970 


3 


-2.925 


-2.936 



Table 1: Threshold chemical potentials for Ag deposition on Au(lOO) in 
absence of surface defects. Different lattice sizes {N x N), and cutoff radii (R) 
for the interactions were employed. 



MCS 


20 


40 


80 


160 


320 


640 


1280 


2560 


Pk 


0.0261 


0.0235 


0.0221 


0.0197 


0.0163 


0.0157 


0.0132 


0.0117 


Ps 


0.108 


0.1044 


0.0991 


0.0919 


0.0877 


0.0716 


0.0694 


0.0560 


p 

A 


82.33 


32.33 


20.27 


12.88 


8.25 


4.12 


3.29 


1.96 



MCS 


5120 


10240 


20480 


40960 


81920 


163840 


327680 


655360 


Pk 


0.0105 


0.0081 


0.0068 


0.0064 


0.0050 


0.0029 


0.0027 


0.0030 


Ps 


0.0432 


0.0392 


0.0353 


0.0260 


0.0184 


0.0172 


0.0145 


0.0142 


P 


1.26 


0.94 


0.77 


0.56 


0.35 


0.27 


0.23 


0.24 



Table 2:Density of kink sites pk, density of border sites at the monoatomic 
steps Ps and perimeter to area relationship ^ of the Au islands generated by 
simulated annealing. MCS is the total number of Monte Carlo steps employed 
for the annealing process between 10.000 K and 300 K. 



MCS 


20 


320 


5.120 


81.920 


655.360 


P^tnfl 


-3.14 


-3.12 


-3.11 


-3.11 


-3.11 


ideat 
Pinfl 


-3.12 


-3.11 


-3.11 


-3.11 


-3.10 



Table 3 Chemical potential pinfi in eV corresponding to the isotherms 
resulting from the simulation of Ag deposition on kink sites. These belong to 
Au imperfect surfaces generated by simulated annealing with different number 
of Monte Carlos steps {MCS). p^^rf are the corresponding values for the 
inflection points of the ideal isotherms. 
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8 Figure Captions 

Figure 1 a) Some of the most common surface defects found on single crystal 
surfaces: 1) flat terrace, 2) monoatomic step, 3) kink site, b), c), d) different 
steps of metal deposition. 

Figure 2 Environment of S sites ( a) S=9, b) S=13) surrounding an adsorp- 
tion node. These were employed for the calculation of the adsorption/desorption 
energies of Ag and Au on Au(lOO) for a cutoff radii corresponding to the distance 
between second (a) and third (b) nearest neighbours. The particle adsorbed is 
located in the central box (9 and 13 respectively) . 

Figure 3 Adsorption energies of a Ag atom at the center of the environment 
shown in figure 2b. The histogram drawn with the continuous line considers only 
the configurations involving Ag atoms, while those configurations involving only 
Au atoms were drawn with the dotted line. The plot drawn with pulses includes 
all the possible configurations that may be obtained with the environment in- 
cluding 13 adsorption nodes. In all cases the proper normalizing factors were 
employed( 4096, 4096 and 531441 respectively). 

Figure 4 a) Adsorption isotherms for Ag deposition on defect free Au(lOO) 
(1-4) and Ag(lOO) (5) single crystal surfaces. Different lattice sizes (N x N) and 
cutoff radii R were employed. In each case R = n denotes that the distance to 
the n*'' nearest neighbour was employed. l)iV — 100, i? = 3; 2) iV = 40, R — 2; 
3) A^ = 40, i? = 3; 4) TV = 100, R ^ 2; 5) N ^ 40 and R ^ 3. 

b) Hysteresis effects for the simulation of Ag adsorption on Au(lOO). The 
simulations were performed on a 100 x 100 lattice. The cutoff radius corre- 
sponded to i? = 3 and the starting conditions for the coverage degree were 
1)6^3 = 1; 2) OAg = 0.5 and 3) SAg = 

Figure 5 Adsorption isotherms for Ag adsorption on Au(lOO) at different 
temperatures. . 

Figure 6 Final configurations for the simulated annealing simulations. The 
number of Monte Carlo steps Nmcs increases from upper left to down right. 
Nmcs = 20 X 2"'~^, where m is the ordinal number of the configuration in the 
figure. 

Figure 7 a) Adsorption isotherms for Ag deposition on kink sites for different 
types of imperfect surfaces obtained through the simulated annealing procedure. 
b)Ideal isotherms calculated according to the procedure indicated in the 
text for the different types of imperfect surfaces. 

Figure 8 Types of kink sites most commonly found after the simulated 
annealing procedure. 

Figure 9 Adsorption isotherms for decoration of monoatomic steps. The Au 
surface structures considered were the final states of simulated annealing runs 
with 20, 320, 5120, 81920 and 655360 MCS per temperature cycle respectively. 

Figure 10 Final state at different chemical potentials for a system in which 
the Au islands were obtained by a simulated annealing process with 5120 MCS 
per cooling cycle. The Au atoms on the surface are represented by black pixels 
and the Ag atoms by grey pixels. 1) /i ~ —3.29 eV, 2) /i = —3.07 eV, 3) 
^l = -2.97 eV, 4) ^i^ -2.94 eV, 5) /i = -2.93 eV, 6) ^l =^ -2.91 eV 
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Figure 11 Adsorption isotherms for Ag deposition on Au(lOO) in the pres- 
ence of surface defects. Each curve corresponds to one of the five systems 
considered previously. 
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